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The hardcore-Bose-Hubbard model with random chemical potential is investigated using quantum 
Monte Carlo simulation. We consider two cases of random distribution of the chemical potential: a 
uniformly random distribution and a correlated distribution. The temperature dependences of the 
superfluid density, the specific heat, and the correlation functions are calculated. If the distribution 
of the randomness is correlated, there exists an intermediate state, which can be thought of as a 
localized condensate state of bosons, between the superfluid state and the normal state. 



Bose-Einstein condensation and superfluidity are two 
important concepts in modern physics. Both of these 
concepts are characteristic of the ground states of quan- 
tum fluids and have been investigated in lots of types 
of systems. However, the relation between Bose-Einstein 
condensation and superfluidity remains unclear [ij. Bose- 
Einstein condensation was flrst predicted by Einstein 
for an ideal Bose gas. The Bose-Einstein condensate 
(BEC) was later observed as the A transition of strongly- 
correlated liquid ^He, and has recently been realized in 
dilute atomic gases [2] and excitons in semiconductors [3]. 
For non-interacting particles, the criterion for Bose- 
Einstein condensation is that the occupation number 
for one of the single-particles energy levels should be 
macroscopic. The general criterion for interacting sys- 
tems is that the one-particle density matrix p(x, y) = 
(■0'l'(x)V'(y)) tends to a constant as |x — y| oo. In con- 
trast, the term superfluidity originally described a group 
of experimental phenomena such as frictionless flow and 
persistent current. After the observation of superfluidity 
in liquid ''He, superfluidity was observed in Fermionic liq- 
uid ^He and electrons in metals as superconductivity, and 
is probably realized even in neutron stars. As proposed 
by Hohenberg and Martin superfluidity is not an equi- 
librium property of the system, but rather an anomaly 
of the transport coefficient. However, our understand- 
ing of superfluidity is considerably less than our under- 
standing of Bose-Einstein condensation. Our first system 
showing Bose-Einstein condensation and superfluidity is 
liquid ^He, in which Bose-Einstein condensation and su- 
perfluidity occur simultaneously at the A temperature, 
but these two concepts are, in principle, not equivalent. 
Recent studies on disordered Bose systems, such as "'He 
in restricted geometry, describe the separation of BEC 
and superfluidity. Furthermore, the nature of both Bose- 
Einstein condensation and superfluidity has been inves- 
tigated. The present study is an attempt to solve the 
essential problem in physics by quantum Monte Carlo 
simulation. 

The disordered Bose system is actualized for liquid 
He confined in the nano-porous glassdiSi and the 
atomic gas in the optical trap potentials [9l|. In particu- 



lar, the confinement of ^He has been investigated exper- 
imentally using variousporous glasses, such as VycorQ, 
Aerogel[10j, and Gelsil[8|. In addition, the confinement 
of the '*He in a uniform-nano-porous medium has also at- 
tracted interest 11 , l2| • The pressure-temperature phase 



diagrams change drastically depending on the dimension- 
alities and pore-sizes of the medium. The superfiuidity 
is suppressed as the pore-size becomes narrow, because 
the interaction between ^He atoms increases with con- 
finement in small pores and applied pressure. For the 
^He in the FSM[11], which is one-dimensional and uni- 
form porous, there exists a critical pressure above which 
the normal-fluid state is maintained and the superfluidity 
does not appear, even at zero temperature, if the diam- 
eter of the pore size is less than 2.8 nm. 

Recently, low-temperature experiments [8^ using ^He in 
nano-porous Gelsil glass, which has a mean pore size of 
2.5 nm, revealed that the onset temperature of the super- 
fluidity was lower than the peak temperature of the spe- 
cific heat that corresponds to the transition temperature 
to the BEC state. This means that between the SF state 
and the normal fluid state, an intermediate state, called 
the localized Bose-Einstein condensate (LBEC) state, ex- 
ists. In the LBEC state, the ^He is in the BEC state, but 
superfluidity cannot be detected through a torsional os- 
cillator experiment. The LBEC state can be thought of 
as a localized state of the BECs by the random potential 
created by the nano-porous medium. Note that for the 
case of the '^He in the uniform medium HMMsflsl, which 
has the same pore size as the mean size of the Gelsil glass, 
there is no clear indication of the existence of the LBEC 
state. In addition, the dependences of the superfluidity 
and the LBEC state on the randomness and pore-size of 
the porous medium are of particular interest. 

The above described experimental results suggest the 
problem of how the superfluidity depends on the ran- 
domness and the finite-size effect of the confinement. In 
addition, "'He in the three-dimensional system of Gel- 
sil glass has different transition temperatures between 
the BEC and the SF state. This can be realized by 
cooperation among the randomness and the finite-size 
effect. This system differs completely from the three- 
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dimensional bulk '*He case, in which the Bose-Einstein 
condensation and SF transition occur simultaneously at 
the A temperature. 

In the present paper, we investigate the randomness 
and its distribution dependence on the Bose system. We 
consider the Bose-Hubbard model with random chemical 
potential. Taking a hard-core limit for simplicity, the 
model Hamiltonian on the three-dimensional cubic lattice 
with periodic boundary condition is given as 



(1) 



where bl(bi) is a boson creation (annihilation) operator, 
fii = blbi, and /i^ is the random potential that depends 
on site i. Here, we define fj,i such that it can takes two 
values, — /Lt- and +^J.+ /i+ > 0). For the random dis- 
tribution of {/ii}, we consider the two cases. One case is a 
uniformly random distribution, that is, the sites at which 
fii = fi^{'-fi-) are uniformly distributed with probability 
P^{P^ = 1 — P+). The other case is a correlated distri- 
bution, that is, the sites at /i^ — — /i_ tend to make clus- 
ters, and the clusters of — /i- arc placed randomly such 
that the density of the sites of — /i- is P_. We set the 
mean size of the clusters to seven to eight sites. We per- 
form quantum Monte Carlo simulation using a directed 
loop algorithm lj| in order to investigate the Hamilto- 
nian ([T]) and calculate the temperature dependences of 
the susceptibility, the SF density, the specific heat, and 
the two-point correlation function. The random average 
for the chemical potential is taken over more than 1,024 
samples. The dilute Bose system 15| and the mean-field 
approach [3| have been investigated in terms of the uni- 
formly distributed random potential at finite tempera- 
ture. The thermodynamic property of the Hamiltonian 
([1]) has not been investigated directly. 

Here, wc introduce two susceptibilities to determine 
the transition temperatures corresponding to the two 
states: the state in which the bosons at /i^ — — /i_ form 
BEC clusters and the SF state over the entire system. 
Note that, according to Harris's criterion 17| , the uni- 
versality class of the three-dimensional XY model is not 
affected by on-site potential. Thus, we can determine 
the critical temperature from the finite-size scaling (FSS) 
using the critical exponents of the three-dimensional XY 
model. We compute the system size dependence of the 
'total' susceptibility as follows: 



,+,(T))rfr, (2) 



where r is the imaginary time. We then perform FSS. 
The critical temperature obtained from the FSS of the xt 
corresponds to the onset temperature of the SF density 
because the correlation length extends over the entire 
system at the temperature. We call this temperature the 
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FIG. 1: (color online). Finite-size-scaling plot of the xi using 
the critical exponents of the three-dimensional XY model, 
7 = 1.3177 and u = 0.67155. The data of the systems of 
different sizes follow the same curves, which means that FSS 
works well. We estimate the transition temperature T\ to the 
locally BEC state in this manner. The error bars of the data 
are smaller than the size of the symbols. 



critical temperature Tt, at which the entire system is in 
the SF state. Next, we define the local BEC order' and 
its susceptibility, as follows: 
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E E f {blm+rir))dr, (3) 



where N_ is the number of sites at which /i = — and 
is a limited summation over these sites. Equa- 
tion ^ is an integration of the correlation between a site 
at which fi = —fi- and a site that is located at a distance 
r distant the fi = —fi- site. From the FSS of the Xh we 
define the critical temperature T\ at which clusters of the 
/ii = — sites are locally in the BEC state. 

Figure [1] shows the FSS plot of xi when //_ = 2.0, 
/i+ = 7.0, and P_ = 0.2, and the random distribution 
has a correlation. Finite-size scaling using the exponents 
of the three-dimensional XY model is found to work well 
for the case in which T\/ J — 0.53. In contrast, the FSS 
of the Xt yields the critical temperature T^/J = 0.49. 
For the case in which the random potential is distributed 
uniformly, there is no difference between T\ and Tt within 
the numerical error. For the case of the correlated ran- 
dom distribution, the bosons at /x; = — locally enter 
the BEC state at T\ and form clusters, which then link 
together as the temperature decreases. When the locally 
correlated clusters percolate through the entire system, 
superfiuidity occurs at Tt. Since we are interested in the 
intermediate state between T\ and Tt, in the remainder of 
the present paper, we discuss only the case in which the 
distribution of the random potential has a correlation. 

Figure [2l^a)shows the temperature dependence of the 
SF density [l9| for various system sizes, and Fig. [Hb) 
shows the SF density as a function of the inverse system 
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FIG. 2: (a) Temperature dependence and (b) the system size 
dependence of the SF density for the case in which \i- = 
2.0, pi+ = 7.0, and P_ = 0.2. The temperatures estimated 
from FSS for Xt and xi are Tt/J = 0.49 and Ti/J = 0.53, 
respectively. The onset temperature of the superfluidity is 
consistent with Tt. The line in (b) is a fitting plot of T/J = 
0.5. (color onUne). 



size XjL. The line in the Fig. [2jb) is a fitting plot of 
T/J — 0.5. The onset temperature of the superfluidity 
corresponds to the critical temperature Tt. Hence we find 
that below Tt both BEC and superfluidity are achieved 
over the entire system. We consider the intermediate 
state of the temperature region between T\ and Tt as 
the state in which the BECs are localized by correlated 
random chemical potential the entire system does not 
show the superfluidity. 

In Fig. [21 we plot the two-point correlation function 



as a function of the distance r. We calculate G(r) for 
system sizes of T = 16 and 20. Comparing the results 
for T = 16 and 20, we find that if r < 6, there is no 
difference between the results. Therefore, in order to 
negate the finite-size effect, we plot only the region of 
r < 6. For the ordinary second-order transition belong- 
ing to the three-dimensional XY universality class, we 
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FIG. 3: Two-point correlation function as a function of dis- 
tance(logarithmic plot). The error bars of the data are smaller 
than the symbols in the figure, (color online). 



have G(r) oc l/ri+'' = l/ri "38 jT|| at the critical tem- 
perature, above (below) which the G(r) decreases expo- 
nentially (converges to some constant value) as r — > oo. 
In the present case, however, the power law behavior is 
observed for a finite temperature range. Because of the 
random potential, the critical behavior is unclear, and 
there appears to be a finite critical region. Since this 
critical region is small, however, FSS analysis is feasible. 
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FIG. 4: Logarithmic plot of the temperature dependence of 
the specific heat. The line in the figure is a fitting plot of 
C (X T^. (color online). 

The temperature dependence of the specific heat is 
shown in Fig. |4l The line is a curve- fitting plot in which 

C - T^. 



The data is well fitted by the line, while for the uni- 
form system /i.^ — fj,, we find the specific heat behaves 
as C ^ T"^ at low temperature as a result of the gapless 
and linear dispersion of the elementary excitations. The 
T-square behavior of the specific heat is also observed in 
the experiment to examine the '''He in the Gelsil glass 13|. 
However, it is thought that T-square behavior in the ex- 
periment is a cross-over of the T-linear behavior, which 
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is characteristic of the glass state. Thus, it is expected 
that C - T as T ^ ofTl. 
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FIG. 5: Phase diagram in the fi- - T plane, (a) An overview 
of the phase boundaries Tbuik, 31, and Tt. The phase bound- 
aries of the system with randomness and non-randomness are 
approximately parallel, (b) An enlarged view of the LBEC re- 
gion. The vertical-striped region is the state where the bosons 
are locally condensed. The error bars of the data are smaller 
than the size of the symbols, (color online). 

Figure [5] shows the phase diagram of the /i_ — T plane 
for fj,+ — 7.0 and P_ = 0.2. In the phase diagram, in 
addition to Tt and Ti, we plot the critical temperature 
^buik of the system with no-randomness. Comparing the 
system with randomness to that without randomness, Tt 
is found to be lower than Tbuiki which means that the 
superfluidity is suppressed because of the chemical po- 
tential Comparing Tt and Ti, Tt is always lower than 
Ti at the same In the vertical-striped region of Fig. 
Eljb), bosons are locally condensed but the entire system 
is not in the SF state. Note that the phase boundaries of 
the system with no-randomness, Tt and Ti, are approxi- 
mately parallel. 

In conclusion, we have investigated a Bose-Hubbard 
model with random chemical potential by quantum 
Monte Carlo simulation. If the distribution of the ran- 
dom potential has correlation, the localization of the 
BECs occurs at a temperature above the onset temper- 



ature of the superfluidity. We have obtained the phase 
diagram of the /i--T plane. In the intermediate state, 
with respect to the distance, the two-point correlation 
function displayed power-law behavior. 

The reader may wonder whether the locally condensed 
state obtained in the present study does not correspond 
exactly to the experimentally observed LBEC stateQ. 
In the present study, the bosons in the cluster of the 
/ii = — /i- first locally enter the BEC state, and the su- 
perfluidity for the entire system then grows as a result 
of quantum fluctuations as the temperature decreases. 
On the other hand, the experiment involving the ^He 
in the various nano-porous media, the critical tempera- 
ture depends on the pore size because of the confinement, 
where the '^He is strongly correlated under the pressure. 
Each cluster has a different critical temperature depend- 
ing on the strength of the correlation. As a result, the SF 
state is experimentally achieved by percolation through 
the fraction of such LBECs. In the present study, an 
intermediate state exists between the locally condensed 
state and the SF state. However, the present situation 
is not equivalent to that of the experiments, because the 
correlation of the boson is too strong owing to the hard- 
core limit. 
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